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Abstract 

This work is devoted to the formulation and derivation of the k — ^/gamma distribution which 
corresponds to A physical fading model. This distribution is composite and is based on the well 
known k — ^ generalized multipath model and the gamma shadowing model. A special case of the 
derived model constitutes the k — Extreme/gamma model which accounts for severe multipath and 
shadowing effects. These models provide accurate characterisation of the simultaneous occurrence of 
multipath fading and shadowing effects. This is achieved thanks to the remarkable flexibility of their 
named parameters which render them capable of providing good fittings to experimental data associated 
with realistic communication scenarios. This is additionally justified by the fact that they include as 
special cases the widely known composite fading models such as Rice/gamma, Nakagami-m/gamma 
and Rayleigh/gamma. Novel analytic expressions are derived for the envelope and power probability 
density function of these distributions which are expressed in a relatively simple algebraic form which is 
convenient to handle both analytically and numerically. As a result, they can be meaningfully utilized in 
the derivation of numerous vital measures in investigations related to the analytic performance evaluation 
of digital communications over composite multipath/shadowing fading channels. 


Index Terms 
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K — /i Distribution, k — n Extreme/gamma distribution, Multipath fading, shadowing, composite fading 
channels, probability. 


I. Introduction 

It is widely known that fading is an effect that significantly degrades communication signals 
during wireless propagation. A common approach for accounting for this phenomenon has been 
viable through exploitation of appropriate statistical distributions. To this effect, statistical models 
such as Rayleigh, Nagakami-m, Weibull and Nakagami-q (Hoyt) have been shown to be capable 
of modelling small-scale fading in Non-Line-Of-Sight (NLOS) communication scenarios whereas 
Nakagami-n (Rice) distribution has been typically utilized in characterizing multipath fading in 
Line-Of-Sight (LOS) communication scenarios, [Hl-lEl and the references therein. Capitalizing 
on the these models, M. D. Yacoub proposed three generalised fading models, namely, the a — ji, 
the n — fi, the rj — fi models and subsequently the X — and the k — n Extreme models, PHl- 
[M. These models have been proved particularly useful due to the remarkable flexibility offered 
by their named parameters which render them capable of providing adequate fittings to results 
obtained by field measurements. Their usefulness is also evident by the fact that they include as 
special cases all the above small-scale fading distributions. 

It is recalled here that a fundamental principle of wireless radio propagation is that multipath 
and shadowing effects occur simultaneously. As a result, in spite of the undoubted usefulness of 
the aforementioned fading models, they all ultimately fail to account concurrently for both shad¬ 
owing and multipath fading. In other words, the characterisation offered by the aforementioned 
fading models is limited to modelling either the one or the other effect. Based on this crucial 
limitation, the need for composite statistical models that can provide adequate characterization 
of the fading effect as a whole, became necessary iQl-llTl. 

Motivated by this, the authors in ifTTl proposed the Rayleigh/gamma fading model, which is 
also known as if-distibution, {K). In the same context, Shankar in [fT^ exploited the flexibility 
of Nakagami-m distribution, which includes Rayleigh distribution as a special case, and intro¬ 
duced the Nakagami-m/gamma composite distribution - or generalised iL-distribution, [Kg]- 
Likewise, the Weibull/gamma composite distribution was proposed in [[T3ll while an introduction 
to generalised composite fading models was reported by the authors in [IT4l - lfT^ . 

The aim of this work is the derivation of novel analytic results for the n — /i/gamma and 
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the n — fx Extreme/gamma composite fading distributions. In more details, after formulating 
these models, novel analytic expressions are derived for their corresponding envelope and power 
probability density function (pdf). The validity of the offered expressions is justified through 
comparisons with results obtained from numerical integrations while their behaviour is assessed 
under different parametric scenarios. Importantly, owing to the relatively simple algebraic repre¬ 
sentation of the offered expressions, they can be considered useful mathematical tools that can 
be efficiently utilized in studies related to the performance analysis of digital communications 
over K — /i/gamma and k — fx Extreme/gamma composite fading channels. To this end, they 
can be exclusively employed in the derivation of analytic expressions for various performance 
measures such as error probability, channel capacity and higher order statistics, among others. 

The remainder of this paper is organised as follows: Section II revisits the basic principles of 
the K — fx, the k — jx Extreme and the gamma distributions. Subsequently, Sections III and IV are 
devoted to the formulation and derivation of the n — /i/gamma and n — fx Extreme/gamma fading 
models, respectively, along with the necessary analysis on their behaviour. Einally, discussions 
on the potential applicability of the derived expressions in wireless communications along with 
closing remarks are given in Section V. 


II. The K-fx, K-fx Extreme and gamma Distributions 
A. The K-fx Fading Model 

The K — fx distribution is a fading model that accounts for the small-scale variations of the 
fading signal in a line-of-sight (EOS) communication scenario and it is written in terms of two 
physical parameters, namely /i and n. In terms of physical interpretation, the former parameter 
is related to the multipath clustering whereas the latter denotes the ratio between the total power 
of the dominant components and the total power of the scattered waves, (H, [[8]|. According to 
m, for a fading signal with envelope R and f = E{R‘^), the k — fx envelope pdf is given by. 


= 


2^{k + 1 ) ^2 


(t)'" Vi (2/rV^(A: + l)^) 


( 1 ) 


where E{.) denotes expectation and f is the root-mean-square (rms) value of R. Eurthermore, 
the In{x) function denotes the modified Bessel function of the first kind and order n whereas 
the terms k > 0 and /i > 0 are obtained by flU 
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EHR^) l + 2k 

^ Var{R:^) (1 + kf ^ ^ 

which is valid for non-negative values. As already mentioned, the K—fx distribution eonstitutes a 
partieularly flexible multipath fading model sinee it ineludes as speeial eases the widely known 
Nakagami-n (Riee), Nakagami-m, Rayleigh and one-side Gaussian distributions, BHl. 


B. The K-fx Extreme Fading Model 

The n — fx Extreme distribution emerges as a notable speeial ease of the n — jx distribution. 
Speeifieally, it is firstly reealled that all the aforementioned multipath models refer exelusively 
to eonditions involving large numbers of paths whieh ultimately meet the prineiple of Central 
Limit Theorem (CLT). Nevertheless, it was shown reeently that the oeeurrenee of large number of 
paths, is not typieally the ease in wireless radio propagation over enelosed environments sueh as 
aeroplanes, trains and buses [HI, [|T9l, [[20ll . The reason underlying this is that sueh environments 
are known to be eharaeterised by only a few number of paths. As a eonsequenee, severe fading 
eonditions, even worse than Rayleigh, are ultimately eonstituted ifT^ . ifT^ . [l20ll . 

By reealling [[7l|, the Nakagami parameter m is regarded as the inverse of the varianee of the 
normalized power of the fading signal and is related to the parameters n and fx by the following 
relationship. 


m = 


fx{l + kY 
l + 2k 


( 3 ) 


At this point, the k — fx Extreme distribution is obtained by keeping m eonstant, assuming 
Kfx ~ 2m and eonsidering a very strong LOS, i.e. fc —)■ oo, and very few multipaths, i.e. /i —?• 0. 
As a result, equation ([B reduees to lITOl eq. ( 16 )], namely. 


Pfl(p) = 


4mJi(4mp) 




g2m(l-|-p2) 

where p = r/f while 5(.) denotes the Dirae delta funetion. 


( 4 ) 
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C. The Gamma Fading Model 

It is recalled that log-normal distribution has been regarded the optimum statistical model 
for characterising the shadowing effect, [[I]|-[|3l|. Nevertheless, in spite of its usefulness, it 
has been largely shown that when it is involved in combination with other elementary and/or 
special function, its algebraic representation becomes intractable. This is particularly the case 
in studies related to the analytical derivation of critical performance measures in the field of 
digital communications over fading channels. Motivated by this, the authors in ifTTI proposed 
the gamma distribution as an accurate substitute to log-normal distribution. Mathematically, the 
envelope pdf of gamma distribution is given by [fTTl eq. (4)], namely. 



(5) 


where the term 6 > 0 is its shaping parameter and The gamma fading model 

has been shown to provide adequate fitting to experimental data that correspond to realistic 
fading conditions. In addition, it is evident that its algebraic representation is quite tractable and 
therefore, easy to handle both analytically and numerically. As a result, it has been undoubtedly 
useful in characterising shadowing and for this reason it has been exploited in the formulation 
of the K and Kq composite multipath/shadowing models, [0111 . IfT^ . 


III. The k - /r/GAMMA Fading Distribution 


A. Model Formulation 

According to the basic principles of statistics, the envelope pdf of a composite statistical 
distribution is deduced by means of superimpose of two or more statistical distributions. In the 
present case, this is realized by superimposing one multipath and one shadowing distribution. 
In mathematical terms, this principle is expressed as. 



( 6 ) 


where I v) denotes the corresponding multipath distribution with mode y. Based on this 

principle, the corresponding k — /i/gamma composite fading distribution is formulated by firstly 
setting r = X and f = y in ^ and substituting in ® along with equation ([5]). To this end, it 
follows immediately that 
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, , 2jj(l + k)‘^xi^ /■°° Vi , 

Py(^) = —:i3T-/ -^-2- dy (7) 

Importantly, the term y‘^ in ([7]) has emerged from the term f ^ whieh denotes that the mean-squared 
value of the fading amplitude is gamma distributed. However, it is noted here that it can be also 
assumed that the root-mean-squared value of the fading amplitude is gamma distributed. In fact, 
this is precisely the difference between the Rayleigh/Lognormal and Suzuki fading models since 
in the former the rms value of the fading amplitude is modelled as being log-normal whereas 
in the latter it is the mean-squared value of this amplitude which is assumed to be log-normal 
Q. As a result, by applying this principle in ([7]), one obtains 


2y{l + kY^^x^ r°° (2/iA/K(l -f K-) ^ j ^ 


B. A Novel Expression for the Envelope pdf 

The derivation of an analytic expression for the envelope pdf of the n — p distribution is 
subject to evaluation of the integral in ([7]). To this end, it is recalled the the Iu{x) function can 
be equivalently represented in terms of the infinite series in ETl eq. (8.445)], and the polynomial 
approximation in |l22l eq. (19)], namely. 


^ ^ r(/ + i)r(i/ + / + i) (2) 

r{n + l) /X 

~ r(/ + i)r(n-/ + i)r(z/ + / + i) U 

Therefore, by making the necessary change of variableill one obtains 


u+21 


straightforwardl)@ 


(9) 

( 10 ) 


4-1 (2/i4«(l + «)^^ - ^ 


r(n -f /)?7,^ 2«^/i+2/ 1^'"2^+^(1 + k)'"2^+^ 

r(/ + i)r(n -1 + i)r(/i + i)y'^+^ 


( 11 ) 


*The polynomial approximation is used in the present analysis. 

^As n —>■ 00, the polynomial approximation reduces to the infinite series. 
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Subsequently, by reealling that r(a;) = (a: — 1)! and substituting in (fTTI) into ([8]), it follows 
immediately that 


Pxi^) = 

1=0 


+ kY+^K^T{n + 

l\T{n - I + l)T{n + l)e^^^n^^-^T{b)n^ 


X 




ye ^dy 


Xi 


( 12 ) 


Importantly, the above integral ean be expressed in elosed-form aeeording to eq. (3.471.9)], 
namely. 


X 


= 2 ( - 1 K,, 


(13) 


where Ky{x) denotes the modified Bessel funetion of the seeond kind, BTI . Therefore, by 
performing the neeessary variable transformation in (fT3l) and substituting in (O yields the 
following explieit expression: 


Px{.x) 



AT^/ . 7\ 1—9/ M+3( + b I 1 1, 

4r(n + /)n^ ^ yU 2 +4 k‘(1 + /s;) 2 +4 

l\T{n - I + l)r(yu + /)e^«r(6)f)^^-5 





/i(l + k) 


+S6{x) 

(14) 


where (5(.) is the Dirae delta funetion and Z = b — y — l + Furthermore, the notation S denotes 
a normalisation sealar eonstant that needs to be determined so that (fT4l) holds as a true pdf. To 
this end, by reealling that p^{x)dx = 1, it immediately follows that 


^ /!r(n-/ + l)r(/r + /)e/^^r(6)0^^-4 io \ V 


//i(l + k) 

7^ 


dx 

(15) 


Notably, the above integral ean be expressed in elosed-form aeeording to Bill eq. (6.561.16)]. 
Based on this, the eonstant S ean be finally expressed as: 


c. ^ r(n + (;, +1) 

^ !!r(t! - 1 + i)ei‘"T(bj2'/‘ ' 

Therefore, by substituting (O into (fT4l) . an explieit expression for the envelope pdf of the k — y 
fading distribution is dedueed, namely. 


















Pxi^) = 1-2^ -L. \ -rTTTTT-X 




^ /!r(n-/ + l)r(6) ''] 21/4(2-1/2 

which to the best of the authors’ knowledge is novel. 


■P/' I A/'l I + b+iJ,+l _ 1 

r(/i + /)(! + K) 2 4(2 2 4 


(17)' 


C. A Novel Expression for the Power pdf 

The corresponding power pdf of the n — p Distribution is straightforwardly dedueed by making 
use of the relationship: Py^iw) = Pp{^/w)/2^/w. To this effect, it follows that 





r{n + l)n^ 

/!r(n-/ +l)r(6) ^ 


^r+hT{p+\) 

21/4(2-1/2 


p-\-Zl-\-b I 1^ p-\-l-\-b _^ 

2/x 2 2 ^Kz 




r(/i + 0(i + «:)-^- 


'4^2 2 


(18) 


Importantly, the algebraic form of the above expressions is relatively eonvenient to handle 
both analytically and numerically. This is particularly useful in deriving analytie expressions for 
various performance measures in teleeommunications where the pdf of fading distributions is 
often involved in complex integrals combined with other elementary and/or special functions. 


IV. The k- p Extreme/gamma Fading Distribution 
A. Model Formulation 

As already mentioned, the n — p Extreme fading model eonstitute a special case of the 
K — p distribution and it aceounts for severe fading eonditions that typieally oeeur in enelosed 
environments. Eikewise the n — p distribution, the n — p Extreme distribution does not account 
for the simultaneous oeeurrence of shadowing. Therefore, by following the same proeedure as in 
the formulation of the n — /i/gamma distribution, the eorresponding k — p Extreme/gamma mul¬ 
tipath/shadowing composite fading model is formulated. To this end, equation (jH) is equivalently 
expressed as: 


Pp{r) = 4me ^4m-j +e ^™5(r) 


(19) 


















Next, by setting r = x and eonsidering that the (rms) of the the fading amplitude is gamma 
distributed, the f term in (fT9l) is represented by y. Henee, by substituting on ® along with dS]), 
it follows that the k — y Extreme/gamma distribution is given by: 


Pxi^) 


e- 2 m 

r{b)Q^ 



Amli 



-\—h 2 m— y. 

y 


dy + 



( 20 ) 


B. A Novel Expression for the Envelope pdf 

The seeond integral in (l20l) is expressed in elosed-form in terms of the gamma funetion. As 
a result, it follows that 


p^{x) = e-^^ + 


4me 


-2m^ -IL j I . X \ 

-2m foo e ye 4m^) 


'VyJ 


T{b)nb 


, 1-6 


dy 


( 21 ) 


X2 


The modified Bessel funetion in can be re-written aeeording to 11221 eq. (19)], namely. 


h 


4m^'] = 




T{n + l)n^ 2 / 2 i-i- 2 i^i-i- 2 «^i+ 2 / 


VyJ ^ r(/ + i)r(n - +i)r(/ + 2 )y^+l 

By substituting the above expression in yields. 


( 22 ) 


Px W 


= e 


— 2 m 




1=0 


m^(^’''i)r(n +/)ni 2i22«-i-3^i-i-2Z 

/!r(n-/ + i)r(/ + 2)r(6)o^ 


y 




y e 


'^dy 


+ S (23) 


X2 


where as in the ease of the k — /i/gamma distribution, the sealar S denotes a normalisation 
eonstant that eonstitutes the derived pdf true. 

Notably, the X 2 integral has the same algebraie form as the Xi integral in (fT^ and therefore 
in (fT3l) . [|^ eq. (3.471.9)]. To this effeet, by making the neeessary ehange of variables and 
substituting in (|2T]) yields the following analytie expression: 


Px{x) 


= 


n 


1=0 


/ 7\ 1_07 ^ b+3l I 15 b+3^ I 7 

T{n + l)n^ 2^2 2 + 4 777. 2 +4 
l\T{n - I + l)r(/ + 2)r(6)fi^+3 



+ s 

( 24 ) 


By reealling that p{x)dx = 0, it follows that 
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^ = 1-e 


—2m 




T'/ I 7 \ 1 —+ — 

T{n + l)n 2 2 + 4777, 2 +4 

^ /!r(n -1 + i)r(/ + 2 )r( 6 )fi^+i 




(25) 


which with the aid of [|2T1 eq. (6.561.16)], yields 

5 = 1 _ e- 2 - - e- 2 - V (& + |) 

^ r(n-/ + i)r(/ + 2)r(6)o-^-3 


(26) 


Finally, by substituting (1261) in (1241) . an analytic expression for the envelope pdf of the k — fi 
Extreme/gamma composite fading model, namely. 


Px{.x) = 1 


T{n + l)n^ 


n 

-2m , 

^^r(n-/ + i)r(/ + 2)r(6) 


X 


^ , 15 
2 2 + 4 K, 


b-i-i 


2x 


Ml*. I 6 

iZ 2 ^^4 


r(6 + y 


2/!iZ 2 ^^4 777, 2 4 




6 -i -3 6-i 3 

2 2 777 2 4 


(27) 


which to the best of the authors’ knowledge, has not been previously reported in the open 
technical literature. 

C. A Novel Expression for the Power pdf 

By recalling that p^{w) = Pp{^/w)/2^/w and substituting accordingly in (l27l) . the corre¬ 
sponding power pdf of the k — mu Extreme/gamma fading model is deduced, namely. 


2m 


, , 1 -e 

PwV^) = ^ ^ 4 


1—21 2m 


r(? 7 ,-|-/)? 7 ,^ 


^ , n 
2 2 ^ 4 K, 


b-l- 


i(2\/t7) Ji¥+lr(6 + y 


r(?T. — /-f i)r(/-f 2 )r( 6 ) ] yjl *’^Y!fz'’ 2*+4777 '’' 2 ^* I 


, - b-i-l 6-1 3 

-^w2 2 777 2 4 

(28) 


V. Numerical Results and Discussions 

This Section is devoted to the demonstration of the general behaviour of the derived analytic 
expressions for the envelope pdf of the n — /r/gamma and k — p Extreme/ gamma fading 
distributions. To this end, Eigure 1 illustrates the pdf of k — /r/gamma with respect to x for 


















II 



Fig. 1. Envelope pdf of the k — p/gamma distribution for b = 1.4, Q = 1.2, n = 2.0 and different values of k 

b = 1.4, = 1.2, jj, = 2.0 and different values of the parameter n. Likewise, Figure 2 eonsiders 
b = 1.4, = 1.2, K = 1.0 and different values of /i. In the same eontext. Figure 3 demonstrates 
the envelope pdf of the k — /i Extreme /gamma distribution for b = 1.2, = 0.8 and different 

values of the Nakagami parameter m. One ean observe the flexibility of the proposed models 
which, as already mentioned, render them capable of providing adequate fittings to experimental 
results. 

A. Usefulness and Applicability in Wireless Communications 

It is widely known that the algebraic representation of crucial performance measures is critical 
in studies related to analytical performance evaluation of digital communications. This is obvious 
by the fact that when the algebraic form of a corresponding measure is convenient, then it is 
ultimately more possible that the derived relationships be expressed in closed-form. Based on 
this, the fact that the form of the offered analytic expressions is relatively simple, renders the 
proposed fading models convenient to handle both analytically and numerically. Therefore, the 
derived expressions can be efficiently applied in various analytic studies relating to the perfor- 
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Fig. 2. Envelope pdf of the k — p/gamma distribution for b = 1.4, Q = 1.2, k = 1.0 and different values of /r 


mance evaluation of digital communications over composite multipath/shadowing fading channels 
including severe fading conditions. Indicatively, they can be meaningfully utilized in deriving 
explicit expressions for important performance measures such as, error probability, probability of 
outage, ergodic capacity, channel capacity under different constraints and higher order statistics. 
It is recalled here that expressions corresponding to the aforementioned measures can obviously 
be derived in both classical and emerging technologies such as single channel and multichannel 
communications, cognitive radio and cooperative systems and optical communications. 

VI. Closing Remarks 

This work was devoted in the introduction, formulation and derivation of the t] — /x/gamma 
and A —/i/gamma fading distributions. These models are formulated from the superimpose of the 
rf — n and A — /x generalised small-scale fading models, respectively, and the gamma shadowing 
model. These distribution are particularly flexible and it that includes as special cases the widely 
known Hoyt, Nakagami-m and Rayleigh fading models. Novel analytic expressions were derived 
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Fig. 3. Envelope pdf of the k — Extreme/gamma distribution for b — 1.2, = 0.8 and different values of m 


for the envelope probability density function which can be considered a useful mathematical 
tool in applications related analytical performance evaluation of digital communications over 
multipath/shadowing channels. 
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